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We discuss the simple microscopic derivation of a hydrophobic effect. Our approach is based on 
the standard functional representation of the partition function of interacting classical particles and 
subsequent passage to collective variables (local densities of the solvent). We get an expression for 
the solvation free energy of solute molecule of any arbitrary shape and derive the nonlinear equation 
for the mean solvent density surrounding the solvated object. We pay a special attention to some 
inconsistencies between the microscopic consideration and the two-length scale mean-field theory 
of hydrophobic effect. 



I. INTRODUCTION 



The key problem in thermodynamics of hydrophobic interactions consists in accounting of the effect of fluctuating 
media on interactions between solvated molecules. In brief, the hydrophobic effect occurs due to expelling the solvent 
from the volume occupied by solute molecule. Hence, the effect of solvation can be accounted by forcing the total 
solvent density to be zero inside the solute molecule. 

Among various theoretical attempts in developing a constructive theory of solvent-solute interactions, special 
attention deserve the works where the hydrophobic effect is treated in the mean-field approximation with explicit 
separation of two characteristic length scales fj, 0, El ■ In this approach the fluctuating solvent density is decomposed 
in two components, the slowly varying field describing the mean solvent (water) density, and the short-ranged density 
fluctuations (usually supposed to be Gaussian) describing correlations in the solvent on scales of order of size of solvent 
molecules. There is a general belief that such an approach is optimal from different points of view: on one hand it 
is physically clear, being " semi- microscopic" , and on other hand, it can be used as a constructive computational tool 
of account of water, much faster than corresponding explicit approaches, but without essential loss of precision. In 
particular, minimization of the free energy in the frameworks of the two-length scale approach allows to determine the 
structure of hydrophobic layer and the solvation free energy of solute molecules of any geometry, if the density-density 
correlation function of the pure solvent is known. 

To be precise, the two-length scale approach deals with the following Hamiltonian 

H [uJ,n\ = ^J w(r)x _1 (r,r')w(r')drdr' + J || (Vn(r)) 2 + W{n(r))} dr + c J n(r)w(r)dr (1) 

where n(r) is the smoothly changing (average) solvent density; w(r) is the field corresponding to the short-ranged 
density fluctuations; x(|r — r'|) is the solvent correlation function in the bulk; a is the phenomenological parameter 
which requires the microscopic determination — see the discussions in ||; the last term describes the interaction 
between short- and long-scale terms with the coupling constant c; and the self-consistent potential W(n(r)) is 
chosen in the common form of the standard Ginzburg-Landau (GL) expansion for the order parameter n(r) as the 
fourth-order polynomial allowing the liquid-vapor phase transition: 

W(n) = ^(n - ni ) 2 (n - n 2 ) 2 (0 < n x < n 2 < 1) (2) 

where ri\ and ri2 are the values of the order parameter n in the vapor and liquid phases correspondingly (below we 
set, if not specified, ri\ = 0) and b is the coupling constant which in combination with the parameter a defines the 
surface tension. 

Following the general scheme of the works 0, we suppose that the solvent cannot penetrate into the volume 
occupied by the solute molecule. Hence the total solvent density p(r) = n(r) + us(r) is nullified inside the solute: 
p(r) = for all r E v in , where Vi n is the volume occupied by the solute molecule. The fact that we force the total 
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solvent density p(r) = n(r) + u(r) to be zero inside the solute, results in an effective interactions between n(r) and 
w(r). In general, one can permit also the direct coupling between n(r) and w(r) everywhere in the solute. Minimizing 
the corresponding free energy functional, we arrive at the set of equations for the profile of the mean solvent density 
n(r) in presence of the solute molecule: 

- aAn(r) + 5Wi ~ n ^ + f dr ' x ^ ( r> r ') n ( r ') = for r e v in (3a) 

SW(n(r)) 

-oAn(r) + g ^ " + U(r) = for r e v out (3b) 

where it is supposed for simplicity that c = (see Q). We are able to compute also the corresponding solvation free 
energy — see [U for details. 

The theory of hydrophobic effect based on the two-scale Hamiltonian with properly adjusted parameters a and 
b describes qualitatively and even quantitatively many physical effects. For example, it reproduces the oscillatory 
behavior in the free energy of interactions of two solute molecules as the function of their mutual distance, and gives 
with good precision the experimentally measured values of solvation free energy of alkane molecules Q . 

However the Hamiltonian is written ad hoc on the basis of physical suppositions on the possibility of separation 
of interactions in "short" and "long" scales. If the two-length scale theory pretends to describe all the peculiarities 
of the hydrophobic effect, it should be confirmed by more solid arguments. For example, one could ask a question 
whether such a theory can be approved microscopically. Just this question is addressed in our paper. The results of 
our attempts to derive the Hamiltonian Q from microscopic consideration are presented in Section IlIII 



II. FROM MICROSCOPIC HAMILTONIAN TO DENSITY FUNCTIONAL 

To make a contents of the paper as self-contained as possible, we describe in this section basic steps of passing 
from microscopic description of classical system with binary interactions to the pre-averaged description in terms of 
collective variables. The hydrophobic effect will be considered in Section InTl 



A. Collective variables 

We begin with a partition function of a grand canonical ensemble for the system of n identical classical particles 

n=0 ' J \i=l J 

The Hamiltonian H n of binary interacting particles has a standard form: 

n n 

ff n =^E/(r,)+X>( r *-^) (5) 

i—l i<j 

where U(r) is an external field and V(r — r') is a pairwise potential. The activity A is related to the chemical potential 
/i by the relation 

. _ 2-KmkT 
' h 2 

and P = l/(kT) is the inverse temperature. 
Let us define the density 



Now we can pass from coordinates of individual particles, (i = 1, n), to the pre-averaged local collective variables, 
p(r) , via the standard technique of constraints in the functional integral 5] . Using the relation 

n 



and the functional Fourier transform 



p(r) -J2 5 ( r ~ r *)l = / V ^ ex P 1 1 / <M r ) P( r ) ~ E 6( - r ~ r *) dr 



(6) 



we get 



oo . . „ / n \ r n 

Z = Y.J V ^^ J lll^ J exp|-t/J$>(rO 



/ 'P{/>(r)[ (-xp<j / i / drp(r)0(r) - / drp(r)J7(r) - | / dr dr' p(r)V(r - r')p(r') }> (7) 



Carrying out the summation over n we obtain 

Z = J J V{p}V{cj)}e 

where 



S(P,0) 



• s '!/'-0} = £ / p(r)V(r - r')p(r')drdr' - /3 A p(r) (V(r) - U(rfj dr - A J e'^^dr 



(8) 



(9) 



(10) 



and <f> rotated to the imaginary axis. After the functional integration over the field p(r), we have 

Z = AT J V{cp}e~~ sw 
where 

S{(f>} = -^ J ^(r)-C/(r))l/- 1 (r-r')(0(r')-C/(r'))rfrrfr'-A J e'^dr (11) 
The equilibrium value of the field <f> can be directly obtained by minimizing the effective action, S, of the system 



8§{ct>} 



(3Xe 



-W(r) 



p / V _1 (r - r') (0(r') - U(r')) dr' = 



The mean density p(r) is determined now by 



1 SlnZ 



/3 J£/(r) 
that reads 

p(r) = J V-\r - r') (0(r') - U(r')) dr' 
After a little algebra we arrive at the nonlinear self-consistent equation for the average density p 

p(r) = A exp |-/3 f[/( r ) + J V(r - r')p(r')dr' 



(12) 



(13) 



(14) 



(15) 



The parameter A and hence, p, is determined by the value of p in the bulk, pb- Rewritten in the terms of bulk density 
(|15|l reads 



p(r) = Ph exp <^ -/3 U(r) + V(r- r') (p(r') - p h ) dr' 



(16) 
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B. Free energy functional 

Equations l|12[l , i|14|) and (|15fl are sufficient to obtain an explicit expression of p and of thermodynamic functionals 
and F. In the saddle point approximation the grand thermodynamic potential f2 is 

" = > Z = T (17) 

Using l(Ti)l and (|T5)l . we obtain 

n = ~\j P( r ) V ( T ~ r')p(r')drdr' -I J p{v)dv (18) 

From (|15|l we can easily get the chemical potential p = i In A + const 

p, = i In p(r) + U(r) + J V(r - r')p(r')dr' + const (19) 

(we can remove the constants in (|19|) . because they determine only the reference state and do not give any contribution 
to physically important quantities). Recalling that 



F = Q + J p{p)p(r)dr (20) 

we arrive at the following expression for the free energy functional, F 

F = \ J P^ V(j ~ r ')p(r')drdr' + ~ J p(r)(lnp(r) - l)dr + J U{v)p{r)dr (21) 

The grand canonical potential f2 determined by (|18|) depends on density distribution p(r) only, while from the 
"thermodynamic viewpoint" it should be the function of p. Taking into account that according to (|19fl the chemical 
potential p in (|18f) is expressed already as a function of p, we can reconstruct the thermodynamically consistent form 
of f2. Starting from the free energy functional l|21|) and applying Legendre transform, we obtain the grand canonical 
potential with explicit dependence on p: 



Q = - J p(r)V{r - r')p(r')drdr' + ^ j p(r)(lnp(r) - l)dr + J U(r)p(r)dr - p J p{v)dv 



(22) 



C. Correlation and response functions 

There are several equivalent ways that allow to obtain the correlation function from density functional 0] . We use 
the procedure based on variation of the density of a liquid with respect to an external potential U(r). According to 
[(| the one-point (direct) correlation function c^^(r) can be defined via the following relation 

p(r) = A exp j-/?[/(r) +c (1) (r)| (23) 
Comparing [12: il) and (|15fl . we arrive at the following explicit expression for c^(r): 

c W(r) = -j3 J V(r - r')p(r')dr' (24) 
Exploiting standard relation (see, for example, 0) 

«"<'■■') <-) 

we get 

c (2) (r,r') = -/3V(r-r') (26) 
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here c' 2 -*(r) is the well-known Orstein-Zernike direct correlation function. The inverse response function %(r, r') 1 
reads 



X(r,r ) 



m-1 5(r-r<) 



c< 2 >(r,r') 



/ \ «(r-r') 



+ pV(r - r') 



(27) 



p(r) p{r) 
The connection to experimentally accessible quantities is established by means of the following equations 

X (|r - r'|) = p h S(r - r') + p 2 h h(\r - r'|); J x (|r - r"!)^^", r')dr" = S(r - r') (28) 

where h(r) = g(r) — 1 is the experimentally measurable correlation function of the pure liquid. 

Let us mention that the results (|27|l 128fl can be obtained by exploiting an exact relation known as the first Yvon 
equation 6] and (|2*8^l . i.e. 



h(r,r') = 



-kT 5p(r) 5{r - r') 



p(r)p(v') SUp) p(r') 



In case of more sophisticated form of entropy term fsee IIII 01 some modifications have to be introduced in 127fl . 
Namely 



, = Po S ( r ~ r ') 

Xlr ' rj p(r)( P0 -p(r)) 

while the equation l|28|) remains without any changes. 



+ 0V(r - r') 



(29) 



III. FROM MICROSCOPIC APPROACH TO GINZBURG-LANDAU-CHANDLER-TYPE THEORY 

In principle, the partition function (0J describes all physical properties of a liquid, including the possible phase 
transitions, while the mean-field free energy functional l|21|l has a single minimum and hence does not describe phase 
transition. This contradiction is due to the point-like nature of particles in our treatment. We can overcome this 
obstacle introducing a finite size for each particle. This is done in the Section IlII CI In the meantime, in Sections 
IIII Al IIII Bl we split the free energy functional (|21|1 in long- and short-ranged parts and derive the solvation free 
energy using the technique of constraints in the functional integral. 



A. Expansion of the free energy 



Consider a solution p(r) = n(r) describing an equilibrium liquid-vapor interface. Supposing that the continuous 
profile n(r) is rippled by fluctuations 0, 0, write the instantaneous value of the density as a sum of two parts: 
p(r) = uj(r) + n(r), where n corresponds to an equilibrium interface and lo - to fluctuations. Expand now the grand 
thermodynamic potential near the profile n(r) in w(r) up to the second order. This is possible if the fluctuations only 
slightly deform the equilibrium profile. Thus, we have 



n\p] = n[n] 



5p(r) 



1 



p ^ {T)dT ' 2 J Sp{r)5p(r') 



s 2 n 



uj(r)u)(r') drdr' 



By definition: 



5Q 



Sp(r) 



0. 



s 2 n 



p—n 



5 P {v)8p{v') 



X -\r-r>) 



p—n 



Substituting (|3"T|) into i|3U|) . we get 



n[p] = Sl[n] + i J w(r)x _1 (r - r>(r')drdr' 



(30) 



(31) 



(32) 
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The density distribution n(r) is long-ranged, so its Fourier image n(k) is nonzero at small k. Using this fact we can 
simplify (|32[1 in the following manner. Consider the part of the Fourier space where n(k) ^ 0. In this region V(k) 
can be expanded in a Taylor series: V(k) ~ v + fcj9,-F(k) + ^kikjdidjV(k) — 1,2,3). Under the symmetry 

arguments the last expression reduces to 



U(k) ~ v Q + ^ 2 k 2 



So, we get 



where 



J n(r)V(r - r')n(r')dr'dr = v J ?i 2 (r)dr + y J 



(Vn(r)) 2 dr 



v =fv(r)dr, v 2 = - J Y 2 V{r)dv (33) 

Thus, fl[n] appears as: 

Q[n] = iy n(r)(lnn(r) - l)dr + y J n 2 (r)dr + y J (Vn(r)) 2 dr (34) 
Collecting all terms we obtain a mean-field Ginzburg-Landau-Chandler-type functional 

V = \J w(r)x _1 (r - r')u{r')drdv' + J {^(Vn(r)) 2 + W{n(r))} dr (35) 



where 

W(n(r)) = —n 2 (r) + in(r)(lnn(r) - 1) - /in(r) (36) 
2 p 

One can clearly see the similarity between l|35f) - ()36|l and CJ-J2J- ^s it has been mentioned in the Introduction, 
the precise form of the potential W(n(r)) is not important - it should only allow the phase transition. Despite 
W(n(r)) does not fit this condition, its simple generalization obtained in Section fill Cl in the frameworks of the same 
formalism satisfies all necessary requirements. The interaction term between fluctuations of short- and long-ranged 
fields identically vanishes in so the coupling constant c in Q should be assigned to zero. 

The difference between the microscopic approach considered in this paper and the approach based on scale separation 
Q becomes much more essential when we put the solute molecule into the solvent described by the Hamiltonian (0) 
and apply the technique of functional constraints as it has been proposed by Li and Kardar and then exploited by 
Chandler et al 0,S||. 



B. The cavitation part of the solvation free energy 



The influence of cavities, extended surfaces and other inhomogeneities of geometric origin immersed in fluctuating 
environment, can be taken into account by means of a technique Q of constraints in functional integration. In brief, 
the presence of a constraint in a partition function generates an artificial " ghost" field coupled with other " real" fields 
describing the fluctuating media. This interaction corresponds to the media (solvent) response on a constraint. In 
particular, the cavitational part of a hydrophobic effect appears due to the requirement to have no solvent particle 
inside the fixed cavity (solute molecule). 

In this Section we start with the microscopic Hamiltonian JSJ and pass to the collective variables (as it has been 
done in section[n| demanding that no particle enter the volume V- ln of a solute molecule. The corresponding partition 
function can be written as follows 



Z= V{p}V{^} J] S[p}e- 3 ^ 



(37) 



where the action S has the form 

(3 



S{ p ^} = il j p(r)V(r-r')p(r')drdr' -fi j p{r) ^(r) - t/(r)) dr - X j e~ p *^dr 

Using again the functional Fourier representation of 5-function and rotation of ip to the imaginary axis, we get 

'P 



Z=Ni j O{0}D{^}exp|| j drdr'(0(r)-C/(r)+r(r)V'(r))F- 1 (r-r')(0(r')-?7(r')+T(r , )V(r')) 



Integrating over the fields p and tp, we arrive at the following expression 



Z = N 2 jT>{(t>}expl^ j drdr / ^(r)-I7(r))v _1 (r-r / )^(r')-^(r / ))+A J dre~ p ^ 

x exp |-| / dvdv' dx" dv'" (4>{t) - U{v^V- x {v - r")U in (r" - r"')V- l {r"' - r') (<£(*)' - tf(r')) 



here 

f dr"M n (r, r")U(r" - r') = S(r' - r) for r, r' G v h 
V in (r,r") = for r G v out 

and integration in l|39|) on r" and r'" actually is carried out in volume ^i n . 
Introducing the new variable 



0(r) = j V(r- r')(<MO - C/(r'))dr' 
the partition function (|39(l can be rewritten as 

Z = M 3 J We-»W 

where 

.§{#} = ^ / 0(r)V(r - r / )0(r / )drtir / - | / 0(r)V in (r,r')i?(r')drdr / 

+ xjcxp^-P \U(r) + J V(t - r') i9(r')dr' 

In a saddle point approximation we have 

=/3 / K(r - r')tf (r')dr' - /3 f ^ n (r, r')#(r')dr' 



- PX J y(r-r')exp|-/3 ^f/(r') + y V(r' - r")tf (r")dr"J j dr' = 
Calculating now the mean density p(r) in a standard manner, we arrive at the following expression 

p(r) = A exp |-/3 ftf(r) + ^ V(r - r')0(r')dr') J 

where $(r) is defined from (|44|) 

0(r) - /" V^r - r')M n (r',r")^(r")rfr" = Aexp |-/3 f tf(r) + Jv(r- r')ti(r')dA j 
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Let us pay attention to the equations l|45|l - (|46(l obtained using the standard technique of functional constraints 
with the microscopic Hamiltonian before separation of the density field p(r) into short- and long-ranged scales. The 
equation ({46(1 is contradictory to some extent. Namely, for all r s Win the left-hand side of (|46|l identically vanishes, 
while the right-hand side is positive. This contradiction can be formally removed by noting that the potential U 
creates an impenetrable wall around the volume v- ln . However, in this case the value of the field $ is undetermined 
inside the volume v- ln - 

So, the account of hydrophobic effect by natural inserting constraint into the functional integral 1|37[) leads to 
internal inconsistency in the two-length scale theory. One can overcome this contradiction by using the two-length 
scale expansion (|35[1 - <|36[) on the first step, and applying the functional constraint on the subsequent step, against the 
direct use of constraint on microscopic level as in 13 7|) . however such a procedure looks a bit artificial. 



C. Solvent particles of finite volume 



The simplest way to take into account the finiteness of the volume occupied by particles of the solvent consists in 
introducing an additional artificial density p v (r) associated with voids between particles, and requiring then the sum 
of two densities, p(r) and p v (r) to be constant everywhere [lOj. Such a condition guarantees that when a particle 
leaves some volume, it is replaced by a "hole". Such consideration is a sort of an off-lattice Flory-Huggins mean- 
field treatment. Following the same steps as in the Section m we arrive finally at the free energy functional, which 
generates the grand potential functional with a double-well structure, describing a phase transition: 

F=\J p{r)V(v - v')p{r')dvdv' + p(r) lnp(r)dr + ~ J (po - p(r)) In (p a - p(r))dr + J U{v)p{v)dv (47) 

where /3 v ( r ) = po — p(r) and po is the inverse volume of a media molecule. Rewriting the entropy term, neglecting 
the meaningless constants and proceeding to the grand potential f2, we get: 

fi = Po I" P^) ln Pi^) + ( 1 _ pM\ m ( ! _ P^l\ dr+ }_ f p( F )V(r-r')p(r')drdr'- [ pp(r)dv+ [ U{r)p{v)dv (48) 
P J Po Po V Po J \ Pa ) 2 J ' J J 



The equilibrium density profile is determined by the condition 



5p(r) 

which leads to the following equation: 



(49) 



p(r) = (p - P(r)) exp j/3 L - U(r) - J V(r - v')p(r')dA J 



(50) 



Let us focus attention on the structure of ft. There are two independent parameters, namely (3 and p, which determine 
the thermodynamic state of the system. In a homogeneous media the density of the grand thermodynamic potential 
may be written as 

P i_ P , . A P\,„ f , P \ , 1 .. -2 



P Po V Po) \ Po) 2 

The quantity vq is determined by Ij33(l . It is easy to see that H51(l describes the phase coexistence when p c = pqVq/2. 
Then the critical temperature T c of the vapor-liquid phase transition is determined by d 2 p/dp 2 — 0, giving T c = 
—vqPq/A. Below T c there are two different phases, while above T c the system is homogenious. Other values of p 
correspond to cases of single phase in the system, vapor or liquid. The value of po is fixed by the bulk density pb- 
Sufficiently far from the solvated object U(r) ~ 0, thus we have: 

In f-M (52) 



„ Po - pb J ^2 

As was mentioned above the parameter po is related to the inverse volume of the molecule. The latter is rather 
arbitrary, because it strongly depends on the convention how the intra-molecular potential is divided into repulsive 
and attractive parts. Let us stress that in our approach we do not use explicitly the pairwise potential of liquid 
molecules V(r — r'), it enters effectively only through the correlation function h(r — r') — see I|27|) -(I28 |) . For example, 
we can associate the parameter po with the repulsive part of potential and the V(r — r') - with the smooth attractive 
part. So, one can say that the relation l|52|) fixes the value of po, while the correlation function of the liquid h(r — r') 
fixes the potential V^(r — r') via l|28|l . 
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IV. DISCUSSION 



The aim of the present paper is two-fold. On one hand, we have derived the self-consistent set of equations which 
determine the solvation free energy of a molecule of arbitrary shape immersed into a solution. The shape of a solute 
molecule is completely defined by the collection of its Van-der-Waals radii and is encoded in the potential U(r) coupled 
to the equilibrium solvent density p(r). In turn, the equilibrium density p(r) is determined by the nonlinear equation 
containing the pairwise density-density correlation function, \( r ~~ r ') of a pure solvent, which is supposed to be 
known and serves as an input in the theory. The full set of corresponding equations l|48l) . JSUJ, (|27[) is collected below 
for convenience: 

F =\J P( v ) v ( v ~ r')p(r')drdr' + ~ J p(r) lnp(r)dr+i J (p - p(r)) In (po - p(r))dr + J £7(r)p(r)dr 



p(r) = (p - p(r)) cxp j/3 J V{v) dr U(r) -Jv(t- r')p(r')dr') } 



(53) 



I; , 1 / /s-i po^(r-r') 



The correlation function %(r — r') and the density pb in the bulk enter as inputs in the theory, while the relation 
between the constants po and pb is established in (|52|l . Let us stress that the free energy of interactions between 
different solute molecules can be straightforwardly computed on the basis of (|53|) . 

On the other hand, we have established the internal inconsistency between the microscopic approach and the 
two-ranged scale Ginzburg-Landau-Chandler (GLC) description of hydrophobic effect. 

Let us remind that the Hamiltonian of the GLC-type theory consists of two phenomenological parts: (i) the 
nonlocal Gaussian term corresponding to short-ranged interactions between pure solvent molecules in the bulk, and 
(ii) the "long-ranged" term coming from smoothly changing profile of "macroscopic" (i.e. averaged over short-ranged 
fluctuations) solvent density taken in a form of a Ginzburg-Landau expansion. These two contributions to the free 
energy are decoupled until the solute molecule is put into a solvent. In the frameworks of GLC theory the solute expels 
the total solvent density p(r) from the volume. As soon as p(r) is a sum of short- and long-ranged contributions, the 
condition p(r) — inside the solvent molecule leads to an effective coupling between short- and long-ranged fields. 
Then one can proceed with the standard thermodynamic formalism and compute the averaged density, the solvation 
free energy, etc. 

The origin of the discrepancy between microscopic approach developed above and GLC-type theory consists in the 
following. If we use the functional constraint nullifying the total solvent density p (as in GLC theory) in the partition 
function with microscopic Hamiltonian (see ea. (|37[l ). then we come to the contradiction: some fields entering in the 
answer cannot be accurately defined (for more details see the last paragraph of Section Till B|) . Only for pure solvent 
the decomposition on short- and long-ranged scales (as in GLC approach) is consistent with microscopic description — 
compare (|35fl - H36[) and Q-©- At the same time the approach described in the current paper does not require any 
coupling between fields describing the microscopic and macroscopic structures of the solvent — in our work we get rid 
of any artificial division of fields in different scale ranges. Moreover, the values of the coefficients a and b in the GLC 
Hamiltonian are not arbitrary, but are dictated by equation l|33|l . 

Meanwhile, we consider it necessary to note that still the two-length scale model, being properly tuned, provides 
rather successful quantitative description of cavitation part of hard spheres and the solvation energies of alkane 
molecules as it has been shown recently, for example, in A]. So, to our point of view the GLC-type approach is 
nevertheless satisfactory as a phenomenological theory. 
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